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We compute exactly the non-equilibrium DC noise in a Luttinger liquid with an impurity and 
an applied voltage. By generalizing Landauer transport theory for Fermi liquids to interacting, 
integrable systems, we relate this noise to the density fluctuations of quasiparticles. We then show 
how to compute these fluctuations using the Bethe ansatz. The non-trivial density correlations from 
the interactions result in a substantial part of the non-equilibrium noise. The final result for the 
noise is a scaling function of the voltage, temperature and impurity coupling. It may eventually be 
observable in tunneling between edges of a fractional quantum Hall effect device. 



I. INTRODUCTION 

The study of fluctuations in systems of non-interacting 
particles is a well-understood problem. Of particular in- 
terest is the noise due to free electrons scattering off an 
impurity, which has been discussed carefully inETH. In 
these calculations, the only interactions are those by in- 
dividual particles with a barrier. The fact that the par- 
ticles do not interact with each other is of course crucial. 
When interactions are present and the system is out of 
equilibrium, it is not possible in general to go beyond 
a perturbative analysis of the fluctuations. In particu- 
lar, the issue of possible singularities in the noise is very 
difflcult to resolvea. 

In this paper, we present an exact calculation of the 
non-equilibrium noise in a one-dimensional (interacting) 
Luttinger liquid with a single impurity. This problem 
has attracted a great deal of attention because the edge 
of a fractionaLquantum Hall device should be a Lut- 
tinger liquidEl^Q. A point contact pinching the sample 
plays the role of the impurity, with the gate voltage al- 
lowing the impurity coupling to be varied. The noise 
here is especially important because DC transport mea- 
surements are not sensitive to the charge of the carriers, 
while non-et[uilibrium noise measurements do depend on 
this chargeQ. We obtain the general current (and voltage) 
noise as a function of the applied voltage, temperature 
and impurity coupling. We check in particular that the 
current noise corresponds to tunneling of uncorrelated 
electrons in the weak-backscattering limit, and to tunnel- 
ing of uncorrelated Laughlin quasiparticles of fractional 
charge in the strong-backscattering limit. 

In a critical theory, conformal invariance usually yields 
the correlators, from which the fluctuations follow. The 
situation changes dramatically when an impurity is intro- 
duced. This adds interactions to the model and destroys 
conformal invariance. One can use the Keldysh formal- 
ism to relate various types of fluctuations to each other, 
and to obtain results to lowest osder in perturbation the- 
ory in the interaction strengthcl. Pejcturbed conformal 
field theory yields other informationQQ (in particular it 
predicts a singularity at finite frequency), but in neither 



approach has it been possible to compute exactly the 
non-equilibrium noise in the presence of interactions. 

In a series of papers with A. Ludwig, we have devel- 
oped a description of a Luttinger liquid in terms of its 
interacting quasiparticles. As in our earlier work, the re- 
sults rely crucially on the f^Qi, that the Luttinger model 
with impurity is integrablelljllj. This formalism makes it 
possible to do exact non-perturbative calculations in the 
presence of the impurity. We have computed the current 
and conductance through tiie impurit300, and the DC 
noise at zero temperaturctS. At zero temperature, the 
fluctuations are pure shot noise coming from the impu- 
rity, so flnding them does not require understanding the 
non-trivial quantum fluctuations of the interacting quasi- 
particles. However, at non-zero temperature, the thermal 
noise and the shot noise can not be separated, so the full 
noise requires understanding these fluctuations. The ef- 
fect of these fluctuations is substantial; without including 
them one does not even obtain the equilibrium Johnson- 
Nyquist formula when the impurity is not present. In 
this paper, we calculate exactly the quasiparticle fluctu- 
ations, and use them to flnd the exact DC noise at any 
temperature, voltage and impurity coupling. 

In sect. 2, we show how to compute density fluctu- 
ations of the quasiparticles in an integrable theory, and 
apply the results to find the current noise in the Luttinger 
model without an impurity. In sect. 3, we couple an im- 
purity to the Luttinger liquid, and find the exact scaling 
curve. This involves a reformulation of Landauer's ap- 
proach suitable for interacting quasiparticles. In sect. 
4, we consider various limits, including the weak and 
strong backscattering limit, and also derive some gen- 
eralized fluctuation-dissipation results. 



II. DENSITY AND CURRENT FLUCTUATIONS 
IN AN INTEGRABLE THEORY 

In this section we first review known results for the 
quasiparticle densities in an integrable theory. We extend 
these calculations to compute the density fiuctuations. 
We then specialize to the case of interest, the Luttinger 
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liquid (in this section, without an impurity). For this 
model, we find the current and charge fluctuations. 



A. General considerations 

We consider an integrable theory whose space of states 
is described entirely in terms of quasiparticles. In an inte- 
grable theory, the scattering is completely elastic and fac- 
torizable, which means that the energy and momentum 
of the quasiparticles are individually conserved, and that 
the multi-particle S matrix decomposes into a product 
of two-particle elements. For simplicity we require that 
the scattering is diagonal as well, which means that the 
types of the particles do not change in the collision. This 
does not mean the scattering is trivial, because there can 
be a phase shift; the S matrix element '5'j^"''^(pi,P2) for 
scattering a particle of type i with one of type j is a non- 
trivial function of the momenta. It will be convenient 
to parametrize the momenta in terms of the rapidity 9 
defined for gapless theories as p oc exp 0. Then scale in- 
variance requires that S matrix element is a function of 
Pi/P2, or equivalently, 9i — 62- (In a massive theory one 
has p oc sinh 6 and relativistic invariance requires also 
that 5 is a function of 61 — 6*2 ; the results of the section 
apply to this case as well.) 

We study the system on a circle of large length L and 
at temperature T. Each quasiparticle has the energy e^; 
in a gapless theory a = ±pi oc exp 6. We introduce 
the level densities ni{0) and the filling fractions fi{9), 
defined so that Lni{0)d9 is the number of allowed states 
for a quasiparticle of type i in the rapidity range between 
and 6 + d9 and fi (9) is the fraction of these which are 
occupied. Thus the density of occupied states per unit 
length Pi{9) is given by 

P{9)=n,{9)M9). 

Requiring that the wave functions be periodic in space 
gives rise to the Bethe equations (setting h = 1, not the 
usual h ^ 1) 
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d9' ^,,{9~9')Pj{9'), 



(1) 



where the kernel is defined by 
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and all rapidity integrals in this paper run from —00 to 
00. These Bethe equations hold for any allowed configu- 
ration, whether it is the equilibrium configuration or not. 
Notice that in a free theory the kernel vanishes and we 
obtain the free-particle relation ni{9) — dpi/d9. 

We generalize the usual analysis by allowing for a 
chemical potential which depends on 9, thus contributing 



exp 



j d9ii,[9)P,{9) 



(2) 



to the Boltzmann weight of a configuration. This gener- 
alization may seem a bit odd, but in an integrable theory 
with diagonal scattering, collisions change neither the en- 
ergy nor type of a particle. Therefore Pi{9) is indepen- 
dent of time, even when a current is flowing. In the next 
section we discuss how the densities change with time, 
but in the very special situation where the change results 
from a single impurity. 

To determine the equilibrium values of the densities, 
we find the configuration which minimizes the free energy. 
This Dsoccdure is known as the thermodynamic Bethe 
ansatia, and is discussed in the appendix. We define 
the "pseudoenergies" by 
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where and the chemical potentials /i^ are scaled by the 
temperature to make them dimensionless. The result is 
the thermodynamic Bethe ansatz equations 



k 



1 



+- 



(4) 



The bar indicates thermal equilibrium values. The equi- 
librium value F of the free energy per unit length is 



F = -T^ / d9 e^{9) In [l + g-^-W+'^'C') 



(5) 



The relations (|) and (|) still hold when the /ii de- 
pends on 9. This makes it possible to compute correla- 
tors by taking functional derivatives of the free energy 
with respect to ^li{9), and then setting Hi to the appro- 
priate value. This is because the partition function can 
be written as the functional integral 



\{'DP^{e) 



-FoL/T+lJ2,J defi,(9)Pi{9) 



where the explicit form of Fq is given in the appendix 
(here, the exact form of the integration measure is some- 
what unclear, but this will not affect the following anal- 
ysis). The equilibrium value —FL/T is the saddle-point 
value of the exponent. Taking derivatives with respect to 
fj-i{9) brings down powers of Pi{9) and gives expectation 
values. For example, the equilibrium particle densities 
are given by 
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ZLSfi,i9) T5n,{9) 
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Of course doing this functional derivation is not necessary 
for determining Pi ; a direct way is to relate it to e using 
the definition (|^) and the Bethe equations (|l|). 
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To find the current noise we will need the density- 
density correlation function, which is given by second- 
order functional derivatives: 
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5^Z 



ZL'^ 6^i,{e)6iJLj{e') 



It is convenient to define the difference of an quantity 
and its equilibrium value, e.g. AP^ = Pi — Pi. We then 
have 



I 



5^F 



LT 5ii^{e)5^ij{e') 



(7) 



In general, this equilibrium value is not diagonal: it can 
be non-zero even when 9^0' and i ^ j. The diag- 
onal correlations come from the Fermi repulsion of the 
particles, while the non-diagonal part arises from the in- 
teractions. However, we show in the appendix that in 
integrable systems there is at least one kind of diagonal 
fluctuation, even if interacting: 



M,{o)Mj{e') ^ 6,,5{e ~ 9') 

It would be interesting to find a deep reason for this 
relation. 

There is a simple alternative formulation of (|^) in 
terms of the q. We discuss this in the appendix . The 
result is written in terms of a function Kij defined by 

K,,{e,9') = $,,(0 - 9') - (1 + e^'('^)-^'('')),5y-,5(0 - 9'). 

(8) 



Then we find that 



APi(0)AP,(0') 



where the inverse is defined via 



fe 



d9" K, 



ik\ 



^S^jS{9-^ 



Unfortunately, this reformulation, while rather elegant, is 
also rather useless because we cannot find a closed-form 
expression for explicitly. The problem in inverting 
is that K is not a function of the difference 9 — 9' but 
of both the variables individually. Observe from that, 
as expected from general considerations, fluctuations of 
densities per unit length are of order 0{1/L). 



B. Current noise in a Luttinger liquid 

We now use these general results to compute the cur- 
rent noise in a particular system, the Luttinger liquid. 
We first find the noise using conformal field theory, and 
then rederive it using the quasiparticle picture. This 
leads to a result which will be crucial to our analysis 



in the next section including an impurity. We find that 
even though the current operator creates multi-particle 
states, only one-particle states contribute to the DC fluc- 
tuations. 

The bulk Hamiltonian of the Luttinger model can be 
written in terms of the charge currentalj 



= -fdx[ji 
9 Jo 



-fn] 



(10) 



The overall constant g follows from the normalization 
of the current; in the fractional quantum Hall problem 
this is fixed by the charge of tha-electron and we have 
g — V when 1/v is an odd integerEl. Since in this system 
the left and right movers are decoupled, we here treat 
only the right movers. We study the model at non-zero 
temperature T, which corresponds to imaginary times on 
a circle of circumference l/T. 

Without an impurity, the charge and current fluctua- 
tions in a Luttinger liquid (or any other conformal field 
theory) can be found very straightforwardly from con- 
formal field theory.Q The left- and right-moving current 
operators jhix -\- 1) and jR{x — t) are individually con- 
served. To find the noise, we need to integrate the two- 
point function of the currents. Conformal invariance fixes 
it to be 



{3B{x,t)jR{S))) = 



[sinh(27r2r(a;-t))]2 



(11) 



and likewise for the left movers (The strange normaliza- 
tion in ( ITTl) arises because of our choice h = 1 instead 
of the more common h — 1. Also, we are ignoring sub- 
tleties coming from the |eljiral anomaly which vanish in 
the infinite- volume limit .liJ) The DC current fiuctuations 
are then given simply by 



{Aiy 



lim 



dt 



:{{jR{^,t),jR{0,0)})+ {L 



(12) 



where we added the contribution of left movers. This is 
the well-known Johnson-Nykqvist formula (for spinless 
electrons), because the conductance G of this Luttinger 
system is equal to g (seeErEl for example). 

We now rederive this result for the current noise from 
the quasiparticles. We have-disfiussed the quasiparticles 
for this system at length intlUla. The values of g where 
the scattering is diagonal are g = l/t, t an integer. We 
therefore restrict to these values, although we expect that 
the final results will extend straightforwardly to general 
coupling. There are t different types of quasiparticles; the 
kink (labeled by -I-), the antikink (— ), and the breathers 
{1 . . .t — 2). Since there is no scale in the theory, these 
right-moving particles are massless, and have the disper- 



sion relation e, 



Pi 



iC , where tua 



m/2 



and rrij — m sin[7rj/(2(t — 1))]. The overall scale m is 
arbitrary and cancels out of any physical quantities. The 
explicit form of the kernels in the Luttinger liquid is 
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discussed for example inli, where we also give a simpler 
form of the thermodynamic Bethe ansatz equations . 

The general fluctuations depend on space and time. To 
compute them using the quasiparticles, one would like to 
be able to recover expressions such as (|l^) based on the 
action of current operators on multiparticle states. Un- 
fortunately the theory, although it appears "almost free" 
from the Bethe-ansatz point of view, is truly interacting, 
and physical observables act on the quasiparticle basis 
in a very complicated fashion. Typically, operators like 
the current or the energy density have, non- vanishing ma- 
trix elements (so-called form factorgl3) between all pairs 
of states with identical number of kinks minus number 
of antikinks; for instance, jfl(x,t) acting on the vacuum 
can create an arbitrary number of breathers and pairs 
of kinks and antikinks. (Operators that are simple in 
the quasiparticle basis, such as the density of kinks, are 
highly non-local in terms of the original physical current 
jn, and presumably of no physical interest.) This makes 
the computation of the space- and time-dependcpt fluc- 
tuations rather difficult, although not impossible.E3 

While local operators are quite complicated, their in- 
tegral over space can become much simpler. A typ- 
ical example is the stress-energy tensor, whose inte- 
gral gives rise to the energy and momentum opera- 
tors, both of which act diagonally on the multiparti- 
cle states. Another example is the C/(l) charge, Q = 
J jndx. With the convention that the original physi- 
cal electrons have Coulomb charge 1, the breathers have 
no charge, but the kink has charge q = +1 and the 
antikink has charge q — —1. In our quasiparticle ba- 
sis, Q acts diagonally, and is a sum of one particle op- 
erators: Q\Oi, . . .9n)q^...q„ = (Z) 9i ) I ^1 ' ' ■ ' ^n)?! ■ ■ ,<Z,i = 

{L J (P+ - P-)d9) \ei, . . . 0n)q,,...,q^. This property holds 
because the sine-Gordon Hamiltonian (and thus its mass- 
less limit) and the charge operator commute, so they can 
therefore be simultaneously diagonalized. It can also be 
painstakingly checked using general famiulas for matrix 
elements of the local current operator.!!!! 

This makes the computation of the DC noise much 
easier than the AC noise. The point is that DC properties 
are obtained by integrating over all times, and since we 
deal with chiral objects (objects which are functions of 
either x — t or x + t), this is the same as integrating 
over space. Thus DC current fluctuations are related to 
charge fluctuations: for instance, as immediately seen 
from (p^), one has (A/)^ = j-(AQ)^ (we set the Fermi 
velocity equal to unity). 

Charge fluctuations can straightforwardly be com- 
puted using the formalism of the previous section. For 
a given configuration characterized by densities Pi, the 
total charge of the right movers is 



Q = L j{P+- P-)d0. 



(13) 



The effect of a finite DC voltage to add to the free energy 
a term —VQ/2. This corresponds to chemical potentials 



of ^+ = = V/{2T) and = for j = 1 . . . t - 

2. Charge fluctuations can now be obtained by taking 
derivatives with respect to V, the rajpidity-independent 
chemical potential. Using (13) with (g) and (0) gives 



Q{v) 



(AQ)2 



-2L 



dF 
W 



d^F do 



(14) 



We can find Q{V) directly from the TBA. By comparing 
the derivative of (^) with (|l|), one finds that when /i is 
independent of 9 



n,{e,V) = T 



86 ■ 



(15) 



Using the definitions of and fi gives 

de±{9,V) 1 



dO P±{e) = T dO 



89 



1 4- e'^±{s,v)-ti± 



The integrand is a total derivative, and from ( p^ we see 
that 



Q{V) = LT 



In 



1 + e-e+-(V/2T) 



We see from |^) that e(oo) — oo, and it is straightforward 
to show that! 



exp [e±(— oo) 



sinh[(i- l)V/{2Tt)] 



smh[V/{2Tt)] 
This gives Q(y) = VL/t and using ( |l^ ) gives 



(Ag)2 = 2LTg 



(16) 



for any V. Using (A/)^ = ■^(AQ)^ reproduces the de- 
sired result (12). 

For completeness, we can use this computation to give 
a quick proof of the diagonal action of Q. For any mul- 
tiparticle state which we denote by \M), characterized 
densities Pi, one has 

(A^l J dxjR{x,t) J dxjR{x,t')\M) = X 

+ J dede'{p+{e) - p^{e)){p+{e') - p^{d')). 

Here X represents the other contributions to the two- 
point function: 



M' 



{M\ / dxjR{x)\M') 



(17) 



and Em is the energy of the multiparticle state \M.)- It 
is convenient to define 
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D{e,e') = A(p+ -p^){e)A{p+ -P-)(6') 

We then have 

{KiY = L j j D{e,e')dede' + 



(18) 



^2gT+^X, 



-X 



(19) 



where the second equahty foUows from the definition of Q 
in (p^, the last equality follows from (^6|), and X has the 
same expression as ( |l7|) but with a sum on intermediate 
states restricted to Em = E'^^j. However, the conformal 
field theory resuh (|l|) says that (A/)2 25T, so X = 0. 
Since each contribution to X is a square modulus, and 
therefore positive, this means that each of these contri- 
butions must vanish individually. Therefore, the compli- 
cated matrix elements of the local current all disappear 
after spatial integration, and the charge operator, which 
is all we need to compute DC current fiuctuations, acts 
diagonally on multiparticle states, and is a simple sum of 
one particle operators. This fact will prove crucial in the 
next section, where we include an impurity. 

We also note that if one were to neglect the non- 
diagonal density fluctuations, one does not obtain the 
Johnson-Nyquist formula even in this simple situation 
with no impurity. This can easily be checked, because if 
we were to treat the particles as free but with a non- Fermi 
distribution function, then K-j^{0,0') = -h5^j5{e ~ 9'). 
Using this to find D{9,9') yields 



A/2 



_ ^^ smhigV/2T) cosh((l - g)V/2T) 
approx-2T ^^^^) • (^O) 



Except at g = 1/2 (where the Luttinger liquid is equiv- 
alent to a free fermion), this contradicts not only the 
Johnson-Nyquist formula but also the simple physical re- 
sult that when there is no impurity, the left movers and 
right movers are not coupled, and the noise should not 
depend on the voltage. 



III. THE EFFECT OF AN IMPURITY ON THE 
NOISE 

In this section we add a single impurity to the Lut- 
tinger liquid. This couples the left and right movers and 
allows backscattejjujig. An impurity at a; = corresponds 
to adding a terrnERO 



Hb = Xcos[(f>L{0) 



to the Hamiltonian; we have written the theory in terms 
of a boson so that Jl = —dx<t>L(x,t) /2'n and jn = 
dx4'R{x,t)/2'K. Basically, the idea is that we express the 
noise in terms of properties of the system far from the 



impurity combined with the S matrix elements for tun- 
neling through the impurity. The former are given in 
the preceding section, while the S matrix elements are 
known.El|The analysis is similar to one given for free 
electrongj; this is possible because of the simple action 
of Q on multiparticle states explained the previous see- 
As previously discussedQ'El, the transport properties 
of a Luttinger liquid with a backscattering impurity can 
be mapped to a Luttinger system of right movers alone. 
The backscattered charge Ql — Qr due to left-right tun- 
neling in the original problem is proportional to the 
charge Q in the new problem. It is not conserved be- 
cause of the impurity scattering. 

To start, we discuss the one-point function of the cur- 
rent in the presfpce of the impurity, quickly rederiving 
earlier resultsOu There are two types of asymptotic 
states; the in states injected by the reservoir from the 
left, and the out states, after scattering through the im- 
purity. The out state must have the same momentum 
and energy because the integrability allows only elastic 
scattering, so a kink can scatter only into a kink or an- 
tikink of the same rapidity. The S matrix element for a 
particle of rapidity 9 and type i to scatter off the impu- 
rity into a particlfi-pf type j is denoted by Sij{9). These 
are written out ir£3; all we will need for the noise are the 
unitarity relation 



k 



kj 



and the magnitude 

T{9) ^ \S++{9)\' 



1 



1 -)_ e2(ff-i)(«-«B)/s 



(21) 



(22) 



where 9b parametrizes the impurity strength and is re- 
lated to \x3 In particular, it follows from the unitarity 
relation that \S+-\^ = 1 - T. 

The effect of the impurity is that the asymptotic states 
1^?)+ and |6')_ are not individual eigenstates of the Hamil- 
tonian any more. The single-particle eigenstates are now 

^ i= {|0)V" + 5++(e)|e)r + 5+_(0)|0)°"*} 
11^)- = ^ {|(?)!" + 5__(0)|0)!"* + S-+{9)\9)1^*} . (23) 

The unitarity of the S matrix yields i{9\\ \\9)j — Sij. 
Describing multiparticle eigenstates is a little more com- 
plicated; for our purpose it will be sufficifiat to think of 
them as tensor products of the states (p3|)E^. 

Although Q acts diagonally on the asymptotic states, 
it acts in a more complicated fashion on the eigenstates of 
the problem with impurity: this time, Q and the Hamil- 
tonian do not commute and cannot be simultaneously 
diagonalized. First one has 



-{9\\Qm+ = l{l + \S++{9)\'-\S+^{9)\' 
= T{9) 



(24) 
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and 



.{e\\Qm- = -ne). 



(25) 



The charge acts on multi-particle states by multiplying 
by the appropriate densities. Since we will be interested 
in the leading contributions to the current, we can neglect 
the effect of the boundary when quantizing the allowed 
rapidities. Therefore the densities are exactly those dis- 
cussed in the previous section. Since the impurity prob- 
lem is again mapped on purely right movers, the relation 
between DC fluctuations of the current and fluctuations 
of the charge holds as before. It follows that 

~I{V) = f [P+ {9, V) - P_ {9, V)] T{9)d9. (26) 



This is the central result oilil'llHl, where it was however 
derived quite differently. 

The current noise follows in a similar fashion. The 
central result of this paper is that 



(A/)2 = L 



D{9,9')T{9)T{9')d9d9' 

[p+(0)(i-7_(f?)) + p_(0)(i-7+(0))] 

xT{9)[l-T{9)]dd, (27) 



where we remind the reader that D{d,d') is the density- 
density correlator defined in (|l^) . Again, fluctuations are 
computed without the impurity, its effect being negligi- 
ble at leading order in L. The first term in ( p7[ ) is the 
equivalent of (^9|), with the T terms arising from ( |2^ ) 
and (pH). The second term is an additional contribution 
at coincident rapidities. It occurs because Q is not di- 
agonal on the eigenstates. Since Q is diagonal on the 
asymptotic states, the only additional contribution due 
to intermediate states is proportional to 

+ {9\\Q\\9)^ ^{9\\Q\\9)+ ^Ti9){l-Ti9)), 

and another with -!-<->—. This correction does not 
occur if the multiparticle state has a double occu- 
pancy at a particular rapidity 9 (a state of the form 
\\ ... ,9,9 ...) -,,,,), because then the intermediate 
state \ \ . . . ,9,9 . . has two particles in the same 
state, which is not allowed by the Bethe equations. This 
gives rise to the additional densities in the second term 
of (E^), whose equilibrium value is 



P+{9)il-f-{9)) + P^i9){l-f+{9)). 

Of course, using the usual arguments (see the appendix) 
the fluctuations in this term are suppressed by a factor of 
1/L, so we can replace it with the product of the individ- 
ual equilibrium values, as in (p7|). Since fluctuations of 
densities per unit length are of order 0(1/ L), both terms 
in ^ are of order 0(1). 

Q acts diagonally on multiparticle states as a sum of 
one particle operators. In the presence of the impurity. 



Q acts again as a sum of one-particle operators, albeit 
non-diagonally. For each rapidity, there are two possible 
states 11^)+ and \\9)-, and in this subspace Q acts as 
a two-by-two matrix. This matrix can alternatively be 
represented by introducing fictitious fermion creation and 
annihilation operators as 



al{9)a+{9) - al{9)a^{9)\ \S+^{9)\ 



(28) 



-al{9)a^{9)Sl^S+- - al{9)a+{9)Sl_S++ 

We can thus represent the DC current as an integral 
over rapidities with (29) as the integrand. This_formally 
coincides with the free-electron calculation in0. There 
are however major differences since populations at dif- 
ferent rapidities are correlated. For completeness, let us 
briefly rederive ( p7[ ) in this approach. Since we neglect 
the (0(1/L)) effect of the impurity on the densities, av- 
erages have to be computed using the asymptotic states 
and the Bethe equations for the bulk system. One then 
checks that for any such states, (a^(0)a_(0)) = and 

{al{9)a+(9l)a-{9')al{9)) = unless 9 = 9'. The first 
term of (E^ follows right away. For the second term, 
we need also to observe that {a\_{9)a+{9)a-{9)a^_{9)) — 
{al{9)a+{9)){a-{9)al{9)) + 0{l/L), and theP+(l-/_) 
and P-{1 — /+) follow for the coincident term. 

We finally discuss how to compute (^7|) numerically. 
The second term is straightforward; the densities follow 
from the TBA equations (^ and the relation (|l5|) . There 
is a trick to make the numerics even easier; as indicated 
in sect. 4.6 below, this piece can be simply written as a 
derivative of the current. The first term is much more 
difficult. Extracting D{9, 9') from (^ is difficult because 
it would require inverting large matrices and taking their 
derivatives numerically. Taking functional derivatives to 
determine D{9, 9') directly from (|^) is also difficult. How- 
ever, if we choose the rapidity-dependent chemical poten- 
tials to be 



M+(0) --M-(e) 



V_ 

2r 



+ x\S++{9)\\ 



this defines a partition function Z(x). Then 



. j j D{9,9')T{9)T{9')d9d9' 



1 d^F 



T dx^ 



(29) 



Using this method, we calculated the noise explicitly at 
g = 1/3. Some curves are given in fig. 1. A more physical 
way to express the impurity coupling is as 6*^ = In [Is /T] , 
where Tb is a crossover parameter analogous to the 
Kondo temperature. (This definition corresponds to set- 
ting the arbitrary scale m = T.) The noise (A/)^ is then 
a function of three variables V, T and Tb. Since (AI)'^ /T 
is diniensionless, it is a function of two variables, say V/T 
and Tb/T. In fig. 1 each curve is a function of Tb /T at a 
fixed value of V/T. We notice that as V/T is4pcreased, a 
peak develops, similarly to the conductances. However, 
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this peak appears at a lower value of V/T than does the 
peak in the conductance. 




In (T_B/T) 



FIG. 1. The noise {AI)^/T as a function of impurity cou- 
pling for four different values of V/T, where the coupling 
g = 1/3. 



IV. VARIOUS LIMITS 

We now consider various limits of ( p7| ) to check our 
expressions and to investigate further physical aspects of 
the noise in the interacting system. A first trivial limit 
is of course the case with no impurity, where T = 1 and 
the results of section 3 reduce to those of section 2. 



A.T = 

At vanishing temperature, there are no fluctuations, so 
D{9, 6') — 0; the noise is pure shot noise. For positive 
voltage, /_ = 0, so the noise reads simply 



(A/)2= p+(0) T{e)-T{e) 



de, 



(30) 



in agreement withEl. There it is also shown that a simple 
fluctuation-dissipation-type relation relates the noise to 
the current. Unfortunately, we have not been able to find 
a finite-temperature analog of this relation. 



B. = 

When the voltage V = 0,we show that the charge fluc- 
tuations D{9, 9') are diagonal and find that the Johnson- 



Nyquist formula holds. We first observe that due to the 
symmetry of the kernels $y = the Bethe equations 
and the TBA equations (H) yield 



at any voltage. Similarly, when I-l+{9) = —fJ.-{9) = lJ-{d), 
we have 

from which it follows that 

n[fi{e)]=n[-^,{e)] e[fi{e)]^e[-fi{e)]. (31) 

We see from section 2 that D{6, 6') can be written as the 
functional derivative 



D{e,e')\ 



L 



5{p+{e)-p^{e)) 



M(e)=o 



(32) 



Because Pi = nf^ and n and e are even in fj,, the only 
non-vanishing term in the functional derivative at /i = 
arises from the direct dependence of the filling fractions 
(^) on /i. Hence, at vanishing voltage. 



Di9,e%^^ = ^nfil-f)5i9-9'). (33) 



Therefore, the charge fiuctuations are diagonal and free 
in the absence of an applied voltage. This fact also follows 
easily from the symmetry of the TBA equations under 
-I- — when V = 0, combined with (S). One then finds 



(A/)2|^^^ = 2 J nf{l - f )Td9 = 2TG, 



(34) 



where we-used the linear-response result for the conduc- 
tance GM Thus in limit F = we recover the Johnson- 
Nyquist formula. This is an extremely non-trivial check 
on our result, because the conductance here depends on 
the impurity coupling and temperature. 



C. Free fermions 

When g = 1/2, the Luttinger liquid is equivalent to 
free fermions. The bulk S matrices are unity, so the 
kernels $y = and the pseudoenergy Tei{9) ~ rriie^ . 



Then K, 
yielding 



9') = -(1/m, 



9') is easy to invert. 



1 



D{9, 9') = - [P+il - /+) + P_(l - /_)] 5(9 - 9') 
and using n+ ~ n_ = n, we have 

y'rT{[7+(i-7_)+7_(i-7+)]T 

- {!+-l.fT^]d9 
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in agreement with the known resultga-S. 

This hmit iUustrates the importance of the fluctua- 
tions. We note that if one were to use this formula for 



f; 7^ 1 /2 as an approximation for the correct formula (27), 
one obtains an answer larger than the correct one, even 
if we use the appropriate interacting distribution func- 
tions derived from the TBA. For example, as T ^ 1, one 
obtains (20) instead of the correct 2gT. The difference 
between the two is substantial for 1/ 7^ 0, especially when 
g is small. 



D. Strong-backscattering limit 



From (|22| ) we see that for 9 « 6b, the transmission 
amplitude T — > 0. Thus the limit Ob 00 is the strong 
backscattering limit, where the impurity splits the sys- 
tem in two. In the expression (|2^) for the noise, in this 
limit we can neglect compared to T, except for very 
energetic particles with 6 k, 6b- However, for 6 ^ 00 
we have e — *■ c» and / — > 0, and the contribution of this 
region to the integral is negligible. We can therefore ap- 
proximate (|2^) in this limit by neglecting the terms 
entirely: 

(a7F« /n[7+(i-7_)+7_(i-7+)] Tdo. 

By simple algebra one finds 

V \ 



- /-) + - /-) = coth ^— j (/+ - /_) 
Hence 

(Alf « coth (^^) (7+ - 7-)^ - coth (^^) 7. 

This is the noise coming from the tunneling of uncorge- 
lated electrons, a result expected on physical grounds!. 



E. Generalized-fluctuation dissipation results 

By the same arguments as in section 2, the conduc- 
tance G = dl/dV can be written as 



(35) 



The voltage difference v between the two channels is 
given bytl 

9 

where Imax is current with no impurity. Since the noise 
in V is proportional to the noise in the reflected (as op- 
posed to transmitted) current, it is therefore given by 



the same arguments of section 3 with the the transmis- 
sion amplitude T replaced with the reflection amplitude 
1 - r. Therefore 



(Av) 



2 — 



L 



D{6, 6'){1 - T{6){1 - T{6'))d6d6' 



+ J [P+{6){l-f_{6)) + P^{6){l-f4e))] 

xT{6){l-T{e))de, (36) 

By simple algebra, and using ( |35| ) , we derive the following 
identity 



(A/)2 = g^Av)^ + 2T(2G - g). 



(37) 



This is the zero-frequency limit of the general relation 
between current and voltage fluctuations, equation (3) 



F. Weak-backscattering limit 

The weak-backscattering limit is given hy 9b —00. 
This is not as easy to treat as the strong backscattering 
limit, because the fluctuations for very-low-energy parti- 
cles are not negligible, preventing a naive expansion in 
1 — T. We therefore have not succeeded in checking the 
weak backscattering limit analytically. Some steps can 
be accomplished however. We can evaluate the behavior 
of the second term in (^ ) using the identity 



T(l -T) = 



9 



d 



2{g - 1) d6B 



T. 



Therefore 



[/+(!-/-) + /+(!-/+)] T{l-T)d6 
f coth(— ) fn(J.-J )J^T 



2(5-1) 

9 . ,V , dl 

2(^^°*^(2T)a^- 



(38) 



In the weak backscattering limit, the current can be ex- 
panded in a power series / = Imax + c{V, g,T)e^^^^^''^" 
(where I^ax = gV), so 



dl 



^ 2{1 - g) {I - Imax) 



and 



(A/)2 « g C0th{V/{2T))iImax - I) + 2T{2G - g) 
D{6, 6'){l - T{e)){l - T{6'))d6d6'. 



On the other hand, in this limit physical arguments say 
that the noise comes from the uncorrelated .tunneling of 
Laughlin quasiparticles of charge g, yielding^ 



(A/)2 « gcoi\i{gV/2T){Imax -I) + 2T{2G - g). 
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This will hold if and only if 

L I I D{e,e'){i-T{e)){i~T{e'))dede' 

V 
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gV 

coth(— -) — coth , „ 



(I. 



I). (39) 



This can be checked explicitly in the case g — 1/2, but 
in general we need to resort to a numerical check. For 
(7 = 1/3 we computed this integral by expressing it as 
the derivative of Z{x), as discussed at the end of sec- 
tion 3. We found Z{x) by numerically solving the TBA 
equations, and took the derivatives numerically. We find 
that (|9|) is indeed satisfied for any voltage as one takes 
9b — > oo. 



V. CONCLUSION 

We have presented in this paper an exact computation 
of the non-equilibrium noise in a Luttinger liquid with 
a single impurity. The approach is somewhat similar to 
the one for free electrons because the model is integrable. 
However, the quasiparticles of this integrable model are 
not free, and the computation requires as a key ingredient 
the equilibrium fluctuations of densities, which are not 
diagonal and are quite complicated. 

The computation of AC properties appears more dif- 
ficult. In addition to the correlations between densities 
at different rapidities induced by the Bethe equations, 
another key ingredient to take into account is that lo- 
cal physical operators are infinite sums of multiparticle 
operators with complicated matrix elements (form fajC^ 
tors). We hope however that a computation similar toc£l 
is possible. 
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APPENDIX A: DENSITY FLUCTUATIONS 
FROM THE THERMODYNAMIC BETHE 
ANSATZ 

The energy per unit length is given simply by summing 
up the individual energies ei{0) of the particles: 



E = J e,{0)P,{9)d9. 



In order to compute the entropy, we utilize the fact 
that solutions of the Bethe equations allow only one 
quasiparticle per level, so that the filling fractions obey 



< /i < 1. The number of states in the interval between 
9 and 9 + d9 is then 

{n,{9)d9)\ 



iPmd9y.i[n,{9)-P,i9)]d9y.. 
Stirling's formula gives the entropy 

S = ^jd9 [n. In 71, - P, InP, - (n, - P,) ln{n, - P,)] . 

i 

A given configuration is completely characterized by the 
knowledge of the densities {Pi{9)}, or by the knowledge 
of the pseudoenergies {ei{9)} defined in (|^). These sets 
are individually complete because the Bethe equation (|l|) 
gives the rii in terms of the Pi, and so knowing the 
determines the Pi and vice versa. Either set can be used 
as convenient variables for the functional integration. In 
what follows we keep dependent variables in intermediate 
computations for notational simplicity. 

The various quantities are obtained by a functional 
integration with the Boltzmann weight e~^^/'^ , where 
F = E — TS. As discussed in sect. 2, the full action 
includes a chemical potential, and it can be written as 



F = 



^ / (e. - /i.) P^d9-TY,J P^ [ln(l 

i 

-He'--^' ln(H-e-''+''')] d9, 



where we generalize the customary situation by allowing 
rapidity-dependent chemical potentials ^i{9). 

The equilibrium values or Pi minimize the action 
( Al]) subject to the constraint of the Bethe equations 
( The first-order variation gives 



-T 



[ln(l 



ln(l 
ln(l + e 



AP. 



Recall that the AP; are given in terms of Ae^ via the 
Bethe equations. Requiring A_F = gives the TBA equa- 
tions (^) with the equilibrium free energy (|^). 

To compute the fluctuations, we expand this action 
to second order around the saddle point values (denoted 
with a bar). Writing A'^^^Pi for the variation of P; in- 
duced by first and second order variations in e we get 

A(2)f/P = ^ j e.A^^^PjT 
/"pe^'-^'ln(l + e-^'+^')Ae^ 
J A(2)p^ [ln(l + e^^-^') + e''-f'' ln(l + e-'^+^'■ 



E 



E 



Pe' 



ln(l 



)] 

(Ae,)2 
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The definition of e relates the fluctuations of n, P and 
e: 



Also, from the Bethe equations (|l|) we have 



Then it follows from the definition of K that 

-5{a - a')6km ^ ~ ( ^ + 'g^ ) ^krnia, a'). 



da 



Using these two identities and the definition of the in- 
verse of K yields the relation (||) relating the density 
fluctuations to K~^. 



It follows that A^^'^Pi depends only on the combination 
Aej + ^(Aej)^. Because the first-order variation AF — 
0, the linear terms in Ae and many others cancel out, 
leaving only 



P, 



(A similajt-result appears in the original paper by Yang 
and Yangl3 on bosons with repulsive delta interactions). 
Hence the fluctuations of e are diagonal: 



{Ae,{e)Ae,{9')) 



1 + e^^'+^' 
LP, ^ 



(Al) 



In terms of the filling fractions fi we can recast this as 

{n0)A,mn,ie')Af,ie')) = ^nj,{l -J,)6,,6{e - 6'). 

This is the same result as one would have for a free 
theory: fiuctuations restricted to the filling fractions are 
free. The difference with a free theory is that we have 
correlated fluctuations (AniArij) as well as (AriiAfj). 

The densi ty-d ensity fluctuation D{6,d') can be ob- 
tained from ( Al) by expressing Ae^ in terms of the AP^. 
The Bethe ansatz equations (|]) yield 



(1 + e^'-^-)AP. + e^'-^'P,Ae, = * 

3 

which can be recast as 

e^"-^'P,Ae, = ^ J K,j{e,a)APjia)dc 



*AP,-, 



where Kij is defined in (g). Plugging this into (Al) gives 
an explicit form of D{9, 9') in terms of an integral of K~^: 



{AP,i9)AP,i9')) 
-tE / daK-\9,a)K-,\9', 



ale 
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